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Abstract 

We investigate the vahdity of the generahzed second law of gravitational thermodynamics 
on the apparent and event horizons in a non-flat FRW universe containing the interacting 
dark energy with dark matter. We show that for the dynamical apparent horizon, the gen- 
eralized second law is always satisfied throughout the history of the universe for any spatial 
curvature and it is independent of the equation of state parameter of the interacting dark 
energy model. Whereas for the cosmological event horizon, the validity of the generalized 
second law depends on the equation of state parameter of the model. 
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1 Introduction 



The present acceleration of the universe expansion has been well established through numerous 
and complementary cosmological observations [1]. A component which is responsible for this 
accelerated expansion usually dubbed "dark energy" (DE). However, the nature of DE is still 
unknown, and people have proposed some candidates to describe it (for review see [2, 3] and 
references therein). 

One of the important questions in cosmology concerns the thermodynamical behavior of the 
accelerated expanding universe driven by DE. It was shown that the Einstein equation can be 
derived from the first law of thermodynamics by assuming the proportionality of entropy and the 
horizon area [4, 5, 6]. In the cosmological context, attempts to disclose the connection between 
Einstein gravity and thermodynamics were carried out. It was shown that the differential form 
of the Priedmann equation in the Friedmann-Robertson- Walker (FRW) universe can be written 
in the form of the first law of thermodynamics on the apparent horizon [7]. Further studies 
on the equivalence between the first law of thermodynamics and Fricdmann equation has been 
investigated in various gravity theories like Gauss-Bonnet, Lovelock and braneworld scenarios 
[7, 8, 9]. 

Besides examining the validity of the thermodynamical interpretation of gravity by expressing 
the gravitational field equations into the first law of thermodynamics in different spacetimes, 
it is also of great interest to investigate the validity of the generalized second law (GSL) of 
thermodynamics in the accelerating universe driven by DE. The GSL of thermodynamics is as 
important as the first law, governing the development of the nature [10, 11, 12, 13, 14, 15, 16, 
17, 18, 19]. 

Here our aim is to investigate the validity of the GSL of gravitational thermodynamics for the 
interacting DE model with DM in a non-flat FRW universe enclosed by the dynamical apparent 
horizon and the cosmological event horizon. Note that in the literature, people usually have 
studied the validity of the GSL for a specific model of DE with special kind of interaction term. 
But we would like to extend it to any DE model with general interaction term. This paper is 
organized as follows. In section 2, we study the DE model in a non-flat FRW universe which is 
in interaction with the DM. In section 3, we investigate the validity of the GSL of gravitational 
thermodynamics for the universe enclosed by the apparent horizon and the event horizon which 
is in thermal equilibrium with the Hawking radiation of the horizon. Also we give an example 
for the case of the cosmological event horizon. In section 4 we investigate the effect of interaction 
term between DE and DM on the GSL. Section 5 is devoted to conclusions. 



2 Interacting DE and DM in FRW cosmology 

The flrst Friedmann equation in FRW cosmology takes the form 

^•2 k Sir , . 

H^ + ^ = ^ {PA + Pm), (1) 



where we take G = 1 and k = 0, 1, —1 represent a flat, closed and open FRW universe, respec- 
tively. Also PA and are the energy density of DE and DM, respectively. 
From (1), we can write 

n^ + nA = i + Qk, (2) 

where we have used the following deflnitions 

_ Snpra _ 87rpA _ k 
i^m- 3^2' ^'A- 3^2' ^''k- a^H''- 
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Since we consider the interaction between DE and DM, pA and pm do not conserve separately. 
Hence the energy conservation equations for DE and DM are 

/9A + 3ff(l+WAW = -Q, (4) 

Pm + S-ffpm = Q, (5) 

where Q stands for the interaction term. For Q > 0, there is an energy transfer from the DE to 
DM. The choice of the interaction between both components was to get a scahng solution to the 
coincidence problem such that the universe approaches a stationary stage in which the ratio of 
DE and DM becomes a constant [20] . The dynamics of interacting DE models with different Q- 
classes have been studied in ample detail by [21]. Here we continue our study without considering 
a specific form for the interaction term. 

Taking the time derivative in both sides of Eq. (1), and using Eqs. (2), (3), (4) and (5), one 
can get the equation of state (EoS) parameter of DE as 

^^ = -3^(^ + ^^ + ^)- 

For the DE dominated universe, i.e. Qa — ^ 1, the above equation yields 

which shows that for the phantom, wa < —1, and quintessence, wa > —1, dominated universe, 
we need to have H > and H < 0, respectively. Note that for the de Sitter universe, i.e. H = 0, 
we have wa = — 1 which behaves like the cosmological constant. 
The deceleration parameter is given by 

« = -(l + f). (8) 
Replacing the term H/H'^ from (6) into (8) yields 

q = Ul + nk + 3nAC0A). (9) 



3 GSL with corrected entropy-area relation 

Here, we study the validity of the GSL of gravitational thermodynamics. According to the GSL, 
entropy of matter and fluids inside the horizon plus the entropy of the horizon do not decrease 
with time [17]. Here like [15], we assume that the temperature of both dark components are 
equal, due to their mutual interaction. Note that Lima and Alcaniz [22] using a very naive 
estimate obtained the present value of the DE temperature as Tj^g ~ 10~^ K. Also Zhou et 
al. [23] estimated the DM temperature as Tj^j^,^ 10~^ K for the present time. Since the 
temperature of the DE at the present time differs from that of the DM, the systems must 
interact for some length of time before they can attain thermal equilibrium. Although in this 
case the local equilibrium hypothesis may no longer hold [23, 24], Karami and Ghaffari [25] 
showed that the contribution of the heat flow between the DE and the DM in the GSL in non- 
equilibrium thermodynamics is very small, 0(10"''). Therefore the equilibrium thermodynamics 
is still preserved. We also limit ourselves to the assumption that the thermal system including 
the DE and DM bounded by the horizon remain in equilibrium so that the temperature of the 
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system must be uniform and the same as the temperature of its boundary. This requires that 
the temperature T of the both DE and DM inside the horizon should be in equihbrium with 
the Hawking temperature associated with the horizon, so we have T = Ti^. This expression 
holds in the local equilibrium hypothesis. If the temperature of the system differs much from 
that of the horizon, there will be spontaneous heat flow between the horizon and the fluid and 
the local equilibrium hypothesis will no longer hold [23, 24, 25]. This is also at variance with 
the FRW geometry. In general, when we consider the thermal equilibrium state of the universe, 
the temperature of the universe is associated with the horizon. 

The entropy of the universe including the DE and DM inside the horizon can be related to 
its energy and pressure in the horizon by Gibb's equation [13] 

TdS = dE + PdV, (10) 

where V = Airlift/ 3 is the volume containing the DE and DM with the radius of the horizon 
and r = Th = l/(27ri?h) is the Hawking temperature of the horizon. Also 

E=^{PA+Pm), (11) 

P = Pf, + P^ = P^ = OJKPK = -—CJAnx. (12) 

Taking the derivative in both sides of (10) with respect to cosmic time t, and using Eqs. (1), 
(2), (3), (4), (5), (11) and (12), we obtain the evolution of the entropy in the universe containing 
the DE and DM as 

S = 37r H^Rl{Ri,Ri, - HRI){1 + Ofe + ^^a^a)- (13) 

Also in addition to the entropy in the universe, there is a geometric entropy on the horizon 
>Sh = vri?^ [13]. The evolution of this horizon entropy is obtained as 

5h = 27ri?hi?h- (14) 
Finally, the GSL due to different contributions of the DE, DM and horizon is obtained as 

5tot = 37rH^RliRi,Ri, - HRl){l + + ^k^k) + ^nR^Rh, (15) 

where ^tot = 5" + S"}! is the total entropy. Note that in Stot we ignored the contribution of 
baryonic matter (BM) ~ 0.04) in comparison with the DM and DE + ~ 0.96). 

Because according to the recent measurements of the supermassive black hole mass function, 
the present entropy of BM S'^jyj = (2.7 ±2.1) x 10^° is five to seven orders of magnitude smaller 
than the DM 8%^ = 6 x 10^^^^ [26]. 

In the next sections, we investigate the validity of the GSL given by Eq. (15) for the 
dynamical apparent and cosmological event horizons. 

3.1 The dynamical appEirent horizon 

The dynamical apparent horizon in the FRW universe is given by [27] 

RA = H-\l + nk)-^/\ (16) 

For k = 0, the apparent horizon is same as the Hubble horizon. 

Recently Cai et al. [27] proofed that the apparent horizon of the FRW universe with any 
spatial curvature has indeed an associated Hawking temperature Ta = \/2'kR\. Cai et al. [27] 
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also showed that the Hawking temperature can be measured by an observer with the Kodoma 
vector inside the apparent horizon. 

If we take the derivative in both sides of (16) with respect to cosmic time t, then we obtain 

. 3(1 + 9., + Qa^a) 

- 2(1 + J^,)3/2 • ^^^> 

Using Eqs. (16) and (17) one can get 

RaRa-HR^- 2i7(i + J],)2 • (18) 

Substituting Eqs. (16), (17) and (18) in (13) and (14) reduce to 

37r 

S = ^^^^ ^ ^^^3 (1 + ^^fc + 3^^ac^a)(1 +^k + ^a^a), (19) 
Stt 

^^ = ^(rTi4F^'+''^+''^"^^- 

Equation (19) shows that for — (^-n^) < ^a < -^(^-nx^)' contribution of entropy of the 
universe inside the dynamical apparent horizon in the GSL is negative, i.e. S < 0. For the late- 
time or the DE dominated universe where J^a ~^ 1 and Ra = H~^, the entropy of the universe 
will be a non-increasing function of time in the quintessence regime with —1 < uj\ < —1/3. 

Equation (20) clears that for wa > — ( '^^^'' ), the contribution of the dynamical apparent 
horizon in the GSL is positive, i.e. Sa > 0. For the DE dominated universe, the entropy of the 
dynamical apparent horizon will be an increasing function of time in the quintessence regime 
with WA > — 1. 

Finally, using Eqs. (19) and (20), the GSL due to different contributions of the DE, DM and 
apparent horizon can be obtained as 

Equation (21) presents that the GSL for the universe containing the interacting DE with DM 
enclosed by the dynamical apparent horizon is always satisfied throughout the history of the 
universe for any spatial curvature and it is independent of the EoS parameter of the interacting 
DE model. 

3.2 The cosmological event horizon 

For the cosmological event horizon defined as 

i?E = o/ - = al (22) 



a 



one can obtain 

Re = HRe - 1. (23) 

For a de Sitter space-time where H is a constant, Eqs. (22) and (23) show that the cosmological 
event horizon radius is and Re = 0. Therefore, in a spatially flat de Sitter universe, the 
event horizon and the apparent horizon, given by Eq. (16), of the universe coincide with each 
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other and there is only one cosmological horizon [11]. For the de Sitter universe, from Eq. (6) 
we have = — 1/J1a hence Eq. (15) shows that Stot = S = 0, which corresponds to a reversible 
adiabatic expansion. 

Substituting Eq. (23) into (13) yields the entropy of the universe inside the cosmological 
event horizon as 

S = -37rij2i2|(l + Ofe + J^awa). (24) 

Equation (24) clears that for a;A < — ( ^q^'° ), the contribution of the cosmological event horizon 

in the GSL is positive, i.e. S > 0. For the late-time universe, the entropy of the universe will 
be an increasing function of time in the phantom regime with ui\ < — 1. 

Following [28], for the quintessence and phantom universe, i?E > and i?E < 0, respectively. 
Therefore, form Eq. (14) one can conclude that for the quintessence universe Se = 2ttReRe > 
and for the phantom universe S'e < 0. 

Substituting Eq. (23) into (15) yields 



'tot 



27ri?E 



RE-^{l + nk + ^AtOA)H^Rl 



(25) 



which shows that the GSL is satisfied, i.e. Stot > 0, when 



2Re /1 + J4\ 
^A<WTPnW7^-{^ • (26) 



3H^R%nA V 

The above constraint on the EoS parameter of the interacting DE model presents that for the 
late-time universe in the quintessence and phantom regime, we have lu\ < ^^^i 1 and 

< - fli - Ij respectively. 

Substituting Eq. (23) into (26) yields uj\ < /{HRe) where 

^ ^^^^^ - -3iP^nZ " v-^T J • 

In Fig. 1, the solid line shows /{HRe) versus HRe for the DE dominated universe where JIa — > 1 
and the dotted line presents /{HRe) = —1/3. Note that for a universe enclosed by the event 
horizon, we have always lo\ < —1/3. Figure 1 clears that for /(HRe) < oja < —1/3 the GSL is 
not satisfied on the cosmological event horizon and it remains valid only for aj\ < /{HRe)- 

Therefore, for the non-fiat FRW universe containing the interacting DE with DM enclosed by 
the cosmological event horizon, the GSL is satisfied for the special range of the EoS parameter 
of the DE model. Contrary to the case of the apparent horizon, the validity of the GSL for the 
cosmological event horizon depends on the EoS parameter of the interacting DE model. 

3.3 A pole-like type phantom universe enclosed by the event horizon 

Here we give an example to study the GSL in the case of the cosmological event horizon. 
Following [28] , we consider a phantom DE model of the universe describing by a pole-like type 
scale factor as 

a{t) = ao{ts - i)"", t<ts, n > 0, (28) 

then one can get 

H = (29) 
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also 



The cosmological event horizon can be obtained as 

dt ts - t 

It 

also 



= = (31) 



Re = 1— < 0. (32) 

n + 1 

Equations (30) and (32) confirm that the model (28) is correspondence to a phantom dominated 
universe. 

The EoS and deceleration parameters of the model (28) are obtained by the help of Eqs. (6) 
and (8), respectively, as 

'^^ = -3^(' + ^^ + ^)' 

q = -l--<-l. (34) 
n 

Equation (33) shows that for the late-time universe, we have a;A = — 1— 1 which is 
the EoS parameter of phantom DE. 

From Eqs. (14), (24), (29), (31), (32) and (33), the entropy of the event horizon and the 
entropy of the universe inside the event horizon can be obtained as 

27r 

Se = 2ttReRe = - , , ^.. {ts -t)<0, (35) 
(n + 1)^ 

^ = r— n3(*^ -*)(-- ^'^y (36) 

Equation (35) shows that the entropy of the event horizon for the model (28) has a negative 
contribution in the GSL throughout the history of the universe. Equation (36) clears that the 
entropy of the universe has a positive contribution in the GSL only when < 1/n. Finally for 
the model (28), the GSL yields 

5tot = .S + 5e = 7)3^ 1 + n^^k) < 0, (37) 

which clears that for the positive spatial curvature, compatible with the present observations 
[29], the GSL breaks down. 



4 The effect of interaction term between DE and DM on the 
GSL 

In our previous analysis, the interaction term between DE and DM is not appeared explicitly 
in the GSL. To see how the DE-DM interaction influences the GSL, wc need to incorporate a 
specific form of the DE model in our analysis. To do this we consider the holographic DE (HDE) 
model which is motivated from the holographic principle [30]. Following [20], the HDE density 
in a closed universe is given by 

PA = Zc^mIL-"^, (38) 
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where c is a positive constant and Mp is the reduced Planck Mass Mp^ = Sir. Recent obser- 
vational data, which have been used to constrain the HDE model, show that for the non-flat 
universe c = O.SlSloiag [31]. Also L is the IR cut-off defined as 

a 

L=-j=smy, (39) 

where y = y/kRE/a. Note that Re is the radial size of the event horizon measured in the r 
direction and L is the radius of the event horizon measured on the sphere of the horizon [20]. 
For the flat universe L = Re- For a specific form of the interaction term between DE and DM 
as Q = 3b^H(p\ + p^) with the coupling constant [32], the EoS parameter for the interacting 
HDE with DM in a non-flat FRW universe is obtained as [33] 

1 2\/0,\ o/l + 
ajj^ = ..-h^i \ 



3 3c V Qa 

where cos y = y^l — c^Q^k/^K- Therefore the coupling constant due to interaction is appeared 
explicitly in the EoS parameter of the HDE. For the dynamical apparent horizon, Eq. (21) 
shows that the GSL is always satisfied throughout the history of the universere regardless of the 
specific form of DE model and interaction term Q. But for the cosmological event horizon, the 
story is different. The GSL for the interacting HDE with DM in a non-flat universe enclosed by 
the event horizon measured from the sphere of the horizon L is obtained from Eq. (15) as 

Stot = SnH^L'^iLL - HL^){1 + J^fe + J^a^a) + 27rLL, (41) 

where the necessary expressions for L and L are given by Eqs. (15) and (16) in [20], respectively. 



as 



Substituting Eqs. (40), (42) and (43) in (41) yields 



L = — (42) 
COS y. (43) 



'tot oTo' 



|2^^(1 + J^ACOs^ y) + {l- ^)0a cosy + 3(6^ -!)(! + Q^) cosyj 



(44) 



which is same as the result given by Eq. (1.6) in [18]. Equation (44) shows that the coupling 
constant 6^ of the interaction term Q does affect the GSL on the radius of the event horizon 
L. For instance, for cosy = 0.99, J7a = 0.73, fife = 0.01 and c = 1 given by [34] for the present 
time, we get 

5tot = ^(6^-0.264), (45) 
which shows that only for 6^ > 0.264 then S'tot > and the GSL is satisfied. 



5 Conclusions 

Here the GSL of gravitational thermodynamics for the interacting DE with DM in a non-flat 
FRW universe is investigated. Some experimental data have implied that our universe is not 
a perfectly flat universe and it possess a small positive curvature [29]. Although it is believed 
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that our universe is fiat, a contribution to the Priedmann equation from spatial curvature is 
still possible if the number of e-foldings is not very large [20] . The boundary of the universe is 
assumed to be enveloped by the dynamical apparent horizon and the cosmological event horizon. 
We assumed that the universe to be in thermal equilibrium with the Hawking temperature on the 
horizon. We found that for the dynamical apparent horizon, the GSL is respected throughout 
the history of the universe for any spatial curvature and it is independent of the EoS parameter 
of the interacting DE model. But for the cosmological event horizon, the GSL is satisfied for 
the special range of the EoS parameter of the model. 

The above results show that the dynamical apparent horizon in comparison with the cosmo- 
logical event horizon, is a good boundary for studying cosmology, since on the apparent horizon 
there is the well known correspondence between the first law of thermodynamics and the Einstein 
equation [35]. In the other words, the Friedmann equations describe local properties of space- 
times and the apparent horizon is determined locally, while the cosmological event horizon, Eq. 
(22) , is determined by global properties of spacetimes [7] . Besides in the dynamic spacetime, the 
horizon thermodynamics is not as simple as that of the static spacetime. The event horizon and 
apparent horizon are in general different surfaces. The definition of thermodynamical quantities 
on the cosmological event horizon in the nonstatic universe are probably ill-defined [17]. 
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Figure 1: /{HRe) — ^^^fii^ ~ 1 versus HRe for the DE dominated universe. The dotted hne 
shows /{HRe) = -1/3. 
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